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While dark matter self-interactions may solve several problems with structure formation, so far 
only the effects of two-body scatterings of dark matter particles have been considered. We show 
that, if a subdominant component of dark matter is charged under an unbroken U{1) gauge group, 
collective dark plasma effects need to be taken into account to understand its dynamics. Plasma 
instabilities can lead to collisionless dark matter shocks in galaxy cluster mergers which might have 
been already observed in the Abell 3827 and 520 clusters. As a concrete model we propose a 
thermally produced dark pair plasma of vectorlike fermions. In this scenario the interacting dark 
matter component is expected to be separated from the stars and the non-interacting dark matter 
halos in cluster collisions. In addition, the missing satellite problem is softened, while constraints 
from all other astrophysical and cosmological observations are avoided. 


INTRODUCTION 

Dark matter (DM) contributes a significant fraction to 
the energy density of the Universe, with an abundance 
about five times as large as that of baryonic matter. How¬ 
ever, our knowledge about this prevalent form of mass in 
the Universe is quite limited. The dark sector could very 
well consist of several distinct species of particles, each 
with their own interactions and dynamics. 

Observations of the Bullet Cluster (IE 0657-558) and 
galactic DM halos suggest that most of DM is collision¬ 
less [1-3]. Nevertheless, there are hints that at least a 
subdominant component of DM might be self-interacting. 
Apart from the missing satellite and cusp vs. core prob¬ 
lems, some observations of galaxy clusters support this 
claim. A recent weak leasing study of the Abell 3827 
cluster [4] discovered an offset between the distribution 
of DM and the visible stars in the central galaxies of the 
cluster, which can be interpreted as a signature of DM 
self-interactions. Similarly, the distribution of DM in the 
Abell 520 cluster has been reconstructed from weak grav¬ 
itational leasing [5, 6], and it has been suggested that a 
mass peak coinciding with the visible hot gas is required. 

It is difficult to conceive how these observations could 
be explained with collisionless DM; and in the case of 
Abell 520 assuming that all of dark matter is self in¬ 
teracting does not result in a correct DM distribution 
either [7]. Instead, an explanation might require that a 
subdominant component of DM behaves similarly to the 
visible ionized gas—developing shock fronts while its ki¬ 
netic energy is converted into heat. In a cluster collision 
this interacting component of DM will then be subject 
to dissipation and consequently slow down, analogously 
to the X-ray emitting hydrogen gas, thus explaining the 
excess of dark mass on top of the visible gas in Abell 520. 

The possible observation of DM plasma is compati¬ 
ble with recently proposed models of a dark disk within 
our own Galaxy [8]. In these models, the galactic dark 
plasma collapses into a thin disk through radiative cool¬ 


ing. In order for this collapse to occur within the lifetime 
of the galaxy, a light “dark electron” is required, which 
is, however, not necessary for generating the mass distri¬ 
bution observed in the Abell 3827 and 520 clusters. 

As a minimal model capturing the essential features of 
plasma dynamics, we study a specific form of DM, a dark 
pair plasma consisting of vector-like fermions and anti¬ 
fermions that interact via dark photons—gauge bosons of 
an unbroken dark C/(l) gauge symmetry. This model can 
explain the existence of starless DM halos in galaxy clus¬ 
ter mergers, as well as solve some well known structure 
formation problems, while it is only weakly constrained 
by Cosmic Microwave Background (CMB) measurements 
and Big Bang Nucleosynthesis (BBN). Effects of energy 
dissipation in the dark sector have been considered pre¬ 
viously [9], specifically in the context of mirror dark mat¬ 
ter [10-12]. Some of the earlier studies of U{1) charged 
DM [13-18] have considered Debye shielding in the dark 
plasma. The possible significance of the Weibel insta¬ 
bility in galaxy collisions was mentioned in [14]. Recent 
advances in the physics of pair plasmas [19-21] allow us 
draw a more dehnitive picture in this work. 

Up to now, studies of DM self-interactions in cluster 
collisions have mainly focused on the effects of individ¬ 
ual scattering events. The aim of the present work is 
to point out that two-body scattering is negligible com¬ 
pared to the collective plasma effects, and that plasma 
physics should be the starting point for understanding 
the phenomenology of charged dark matter. We argue 
that models where all of DM is charged are ruled out 
by observations of cluster collisions unless the charge is 
extremely small. By keeping this in mind we propose a 
scenario where only a subcomponent of DM is charged. 

A dark pair plasma is an appealing form of DM both 
theoretically as well as experimentally. We show that 
it can be a thermal relic, avoiding the complicated pro¬ 
duction mechanisms required for asymmetric DM such 
as atomic DM [22-26]. It does not generate dark acous¬ 
tic oscillations (DAO) at observable scales, is consistent 
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with CMB measurements, and its dark halos do not col¬ 
lapse to disks, due to inefficient dark bremsstrahlung. In 
addition, dark U{1) models may solve one of the major 
unexplained puzzles of the standard model (SM)—the hi¬ 
erarchy of Yukawa couplings that spread over at least 6 
orders of magnitude [27]. 


DARK PLASMA DYNAMICS 
A Minimal Model 

Plasma is a state of matter in which collective effects, 
mediated by long range interactions, dominate over hard, 
short-range collisions of particles. The minimal model of 
a dark plasma is a pair plasma consisting of thermally 
produced vector-like dark fermions and anti-fermions. 
The Lagrangian of the dark sector is dark QED: 

C = CsM - ^FDt_c„F^‘' + x{0- mo) X, ( 1 ) 

where is the covariant derivative, 

F^'' is the field tensor of the dark photon y is the 
interacting DM component with mass mo and eo is the 
dark U{1) charge. 

This Lagrangian can be extended by including a ki¬ 
netic mixing term Fu^^F^^ that results in electrically 
charged DM. Such a term is severely constrained by re¬ 
combination and halo dynamics [15]. Therefore, based 
on experimental results, we assume it to be negligible for 
the rest of this work.^ 

Nevertheless, we require an interaction between the 
two sectors at some high scale to bring them into ther¬ 
mal equilibrium. It is in principle enough to just assume 
that the visible and hidden sectors are connected by their 
coupling to the inflaton, and therefore share a common 
reheating temperature. Here we will omit the details of 
this interaction, since they are irrelevant to the purpose 
of this paper. Instead we refer to [ 8 ] where various pos¬ 
sibilities for generating the coupling between the visible 
and hidden sectors without inducing the kinetic mixing 
term are reviewed. 

We then assume that the dark sector is populated in 
the early Universe through freeze-out of some feeble inter¬ 
action with the SM above the electroweak scale. There¬ 
fore the temperatures of the two sectors coincide at this 
time, but will evolve differently as the relativistic degrees 
of freedom drop out of equilibrium in the two sectors.The 


ratio of the dark sector temperature To to the photon 
temperature is fixed by entropy conservation 

g*sn(Ty)/9*s,jiT*) 
\g*s,D{TD)/g>,s,D{T*)) 

where and g*sp are the numbers of relativistic de¬ 
grees of freedom in the two sectors, and T* is the tem¬ 
perature at which the dark and visible sectors were pre¬ 
sumably in thermal equilibrium. 

The relic density of the fermionic interacting DM is 
fixed through freeze-out of the annihilation into dark 
photons. The thermally averaged cross section for the 
process XX IdId in the limit Tp ^ mp is 
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where ap = e|,/ 47 r is the fine structure constant of the 
dark U{1). The Sommerfeld enhancement of this cross 
section has a negligible impact on the abundance of the 
dark fermions [13] and will be ignored. By solving the 
Boltzmann equation (using the procedure described e.g. 
in [28]) in terms of the dark sector temperature and ex¬ 
pressing the final result in terms of the temperature of 
the visible sector we can estimate the relic abundance of 
the interacting DM component as 


n 3 ^ ( f ™ V ( ( 4 ) 
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where apM = 1/137, x = mp/Tp and 5 * is the effective 
number of relativistic degrees of freedom in both sectors, 
evaluated at the time of freeze out. The freeze-out tem¬ 
perature can be approximated by 
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and the ratio of the hidden and visible sector tem¬ 
peratures at freeze out assumes values in the range 
Cf € (0.5,1.5), depending on the fermion mass. 

We set this relic abundance to = /Dcdm, 
where ^ is the fraction of the interacting species and 
h^XtcuM = 0.1198 ± 0.0015 [29] is the overall DM abun¬ 
dance. Thus, for 5 = 30% interacting DM, the expected 
relic abundance is obtained approximately if 
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^ To avoid generating the problematic kinetic mixing term radia- 
tively [ 16 ], we have to assume that particles charged under both 
the hidden C/(l) and the SM hypercharge do not exist, or that 
they only exist in complete non-degenerate multiplets of a uni¬ 
fying gauge group. In this case, setting the tree-level kinetic 
mixing term to zero does not require any finetuning. 


Collisionless Shocks 

Collisionless shocks are a prevalent phenomenon in as- 
trophysical plasmas [30], and have been observed e.g. in 
the Earth’s bow shock, in the expansion of supernova 
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remnants into the interstellar medium, and in the behav¬ 
ior of the ionized gas in galaxy collisions. In all those 
situations, the mean free path of particles in the plasmas 
is orders of magnitude larger than the physical size of the 
shock fronts. 

The physics of collisionless shocks, which arise from 
collective plasma instabilities, is an active research topic. 
Laboratory experiments and computer simulations have 
investigated the formation of instabilities in relativistic 
and non-relativistic plasmas consisting of electrons and 
protons, and of electron/positron pairs. 

The formation of collisionless shocks can be roughly 
divided into two phases [20]. The first phase consists 
of the the buildup of the instabilities, followed by their 
saturation. The buildup phase can be studied by linear 
approximations and is relatively well understood. De¬ 
pending on the type of plasma different instability modes 
can dominate [19]. 

In the counter-streaming situation that is relevant 
in cluster mergers, initial fluctuations in the counter¬ 
streaming electric currents in the plasma give rise to 
magnetic fields. These fields enhance the electric cur¬ 
rents, which in turn enhance the magnetic fields, leading 
to an exponential growth. Finally the field strength will 
be large enough to stop the stream of incoming particles, 
and a shock front develops. 

The second phase, where the electromagnetic fields and 
currents reach saturation, is highly nonlinear and can be 
studied only by means of numerical simulations. Most of 
the dissipation of kinetic energy takes place in this phase. 
The latter can be roughly understood as an effect caused 
by scattering of the upstream particles from the strong 
electric and/or magnetic fields generated during the first 
phase. 

The effect of the collisionless shocks in a cluster merger 
can be understood as follows. The DM halos of the 
colliding subclusters are initially in a stable equilibrium 
state. Once they begin to overlap, an unstable counter¬ 
streaming situation arises, and a shock front quickly de¬ 
velops as described above. This shock front then prop¬ 
agates through the subcluster halo, heating and slowing 
down the interacting component of DM, similarly to what 
is seen in the X-ray emitting visible gas. Consequently, 
a generic prediction of the dark plasma model is the ex¬ 
istence of the bow-shaped dark matter shock fronts, that 
should be visible in the weak leasing reconstructions of 
the cluster merger events, given sufficient angular reso¬ 
lution. 

Our suggestion is that these effects might have already 
been discovered in the Abell 520 and 3827 clusters. In 
Abell 520 an excess of DM on top of the visible X- 
ray emitting plasma between the subclusters is observed. 
This can be interpreted as the interacting component of 
dark matter that was slowed down due to the shocks, 
similarly to the visible plasma. 

In Abell 3827 the separation between the stars and the 


center of mass of the dark matter in the central galaxies 
can be interpreted along the lines of [31]: The interact¬ 
ing component of DM in the galaxies is counter-streaming 
against the DM in the main cluster halo. The resulting 
shocks create an effective drag force slowing down the 
interacting DM component of the galaxies, resulting in 
a separation between this component and the rest of the 
mass. Therefore, the center of mass of the total dark mat¬ 
ter distribution is separated from the stars, even though 
the main component of DM remains on top of the stars in 
this scenario. Given a high enough resolution, this effect 
should be observable as separated starless dark matter 
clumps, similar to what is observed in a larger scale in 
Abell 520. 


Shock Formation Time Scale 


To show that these effects indeed are relevant in a typi¬ 
cal cluster merger, we shall now examine the fundamental 
characteristics of the dark plasma—its Debye length, the 
plasma parameter and the plasma frequency. For the col¬ 
liding intracluster DM halos, we set the size to R = 200 
kpc and the mass to M = 4 • 10^^Mq, corresponding to 
the dimensions of the colliding Abell 520 subclusters as 
determined from weak leasing data [6]. Assuming uni¬ 
form distribution, the average density of the interacting 
DM is then 1.36 • 10“^ GeV/cm^. 

This information, together with the interaction 
strength a^i, is sufficient to calculate the time scale for 
the plasma instabilities. The mean free path of DM par¬ 
ticles depends on the temperature of the self-interacting 
DM component, which we can estimate from the virial 
theorem for each of the colliding sub-clusters 
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Here is Newton’s constant, and Udof = 3 is 

the number of degrees of freedom of a single particle. 

Returning to the dark plasma, the Debye length in our 
minimal model is 


If two particles are separated by more than this dis¬ 
tance, the Goulomb interaction between them is effec¬ 
tively shielded by free charge carriers. Clearly this dis¬ 
tance is diminutive in comparison with astrophysical 
scale. 

The plasma parameter A is defined as the number of 
charge carriers within a sphere of radius Ad, 
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Since A ^ 1, the the plasma is weakly coupled and collec¬ 
tive effects caused by long range forces dominate plasma 
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dynamics. The characteristic time scale for collective 
plasma effects is the inverse of plasma frequency 


1/Wp = 


/ mp \V2 
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( 10 ) 


where n is the number density and niD the mass of the 
DM particles. 

The mean free path of charged particles in plasma is 
of the order of [32] 
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The time it takes to form a collisionless shockwave 
can be estimated by considering the instability growth 
rate. In a symmetric non-relativistic collision of two cold 
counter-streaming pair plasmas the dominant instability 
mode is the Two-Stream mode for which the instability 
growth rate is of the order of the plasma frequency [19]. 
A realistic description of the collision certainly requires 
a more complicated set-up than a simple cold counter¬ 
streaming plasma as the latter does not account for the 
possible inhomogeneities and temperature or velocity dis¬ 
persion inside the cluster. However, the latter consider¬ 
ations should not significantly affect the order of magni¬ 
tude estimate given here. A conservative order of mag¬ 
nitude estimate of the time scale of shockwave formation 
is then 


T, « « 57.2 s (12) 

The distance for which plasma instabilities become rel¬ 
evant can be estimated by multiplying our estimate for 
the shock formation time in eq. (12) by the typical speed 
of a cluster collision, 

As « TsUcoi ^ 10^ km, (13) 

which is certainly much smaller than the mean free path. 
Thus, the plasma instabilities affect the dynamics of the 
interacting DM component at time and distance scales 
much shorter than the two-body scattering processes. 
This leads us to conclude that DM charged under an un¬ 
broken [/(I) interaction should be treated primarily as 
a plasma that develops collisionless shocks at relatively 
small scales. Consequently, at larger scales it behaves ef¬ 
fectively as a collisional fluid, even if the two-body scat¬ 
tering rate seems to suggest that collisions are insignifi¬ 
cant. 

We are not aware of studies of non-relativistic pair 
plasmas. However, a study of relativistic pair plasma col¬ 
lisions [21] suggests a time scale that is even an order of 
magnitude smaller than estimated in eq. (12). Our argu¬ 
ment here is based on the instability growth rates in the 
theoretically well understood linear regime. A numerical 
study, which will provide a more conclusive treatment, is 
under preparation. 


Atomic Dark Plasma 

Our minimal model assumes that the interacting sub¬ 
component of DM exists in the form of a collisionless 
pair plasma. It is also conceivable that the dark plasma 
could consist of light “dark electrons” and more massive 
“dark protons”, thus imitating the visible sector. Such 
scenarios of atomic DM require a non-thermal, asymmet¬ 
ric production history. Here we do not discuss any such 
models in detail, but note that any mass non-degeneracy 
of the particles making up the plasma does not signifi¬ 
cantly affect the estimate of the instability growth rate 
[19]. Therefore the collisionless shock behavior will also 
dominate the dynamics of atomic dark plasmas. 

An interesting possibility for explaining both the sub¬ 
dominant interacting DM species as well as the collision¬ 
less main component could be a model of partially ionized 
atomic dark matter, where the dark plasma consists of 
dark ions and dark electrons, while the main component 
of DM is composed of the neutral atoms. 


OBSERVATIONAL CONSTRAINTS 
Bullet Cluster 

This work was motivated by observations of the Abell 
3827 and 520 clusters, but many similar systems exist 
[33]. Perhaps the most unambiguous of those is the Bullet 
Cluster, from which constraints have been derived on the 
interactions of DM. It has been shown that no more than 
30% of the total DM mass can be lost from the subcluster 
as it passes through the main halo [1]. Thus, throughout 
this work we will assume that the fraction of interacting 
dark matter species obeys ^ < 0.3. This rules out models 
where all of dark matter is charged under an unbroken 
1 /(1) as long as the coupling constant is not negligibly 
small. 


BBN 

Changing the properties of DM changes the dynamics 
of the early Universe, and therefore any deviations from 
the usual ACDM model will be strongly constrained. 
New relativistic degrees of freedom change the expan¬ 
sion history of the Universe during BBN, leading to a 
constraint which is usually expressed as a limit on the 
effective number of light neutrino species [29], 

Neff = 3.04 + 2Cdbw = 3.15 ±0.23, (14) 

Assuming one species of fermions we see that the ratio 
(2) of the temperatures of the dark and visible sectors at 
BBN is Cbbat = 0.52, implying that = 3.18. There¬ 
fore Eq. (14) is satisfied within Icr. 
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By considering an extended scenario with No effective 
relativistic dark fermions at T*, the bound in eq. (14) 
implies 


iVD = 0.68 ± 1.67. (15) 

Note that if all fermion masses exceed T^,, then Njj = 0 
and eq. (14) is always satisfied. 


CMB 


Dark matter is usually considered to be pressureless, 
so that all primordial DM density fluctuations start to 
grow immediately after they enter the horizon. If a sub¬ 
component of DM interacts with massless dark photons, 
the growth of structure is suppressed until the DM and 
dark photons kinetically decouple. This effect could be 
observed in the Cosmic Microwave Background as a sup¬ 
pression of fluctuations at large multipole moments. 

The kinetic decoupling of the dark fermions and the 
dark photons occurs when the Compton scattering rate 
in the dark plasma drops below the Hubble rate [34]. The 
Compton scattering rate for the dark plasma is 


^ 135m|) 


(16) 


and the Hubble rate in the radiation dominated epoch is 


H = 




mp ’ 


(17) 


where g* is the effective number of relativistic degrees of 
freedom in the visible sector, with g* = 3.36 at tempera¬ 
tures well below the electron mass, and we have neglected 
the subdominant contribution of the colder dark sector. 
Setting Tc = H we obtain the temperature of kinetic 
decoupling 


on the size of the gravitationally bounded DM struc¬ 
tures due to kinetic coupling with the dark radiation 
is given by ~ 10“^(Tkin/10 MeV)“^M 0 [36], so that for 
Tkin « 0.5 keV the cut off is at ~ 10^Mq, as required to 
ease the missing satellites problem. However, in our case 
only a subdominant part of DM is coupled to the dark 
radiation and thus the cut-off will only affect the inter¬ 
acting fraction of DM, so that structures smaller than 
the cut-off will still exist, only in fewer numbers. 

Figure 1 depicts the kinetic decoupling constraint on 
the (mu, Q;D)-plane, as well as the parameter space re¬ 
gion compatible with the relic abundance. Restricting to 
the region that produces the desired relic abundance, the 
kinetic decoupling limit gives a lower limit on the mass 
of the DM particle. For ^ = 0.3 this limit is roughly 
mp 5 MeV. We show also the upper limit from requir¬ 
ing that the Landau pole of the dark U{1) lies above the 
Planck scale, which for ^ = 0.3 corresponds to roughly 
mp ^ 2 TeV. 

It should be stressed that the constraints for dark mat¬ 
ter self interactions, including the kinetic decoupling con¬ 
straint depicted in the figure, are in reality functions of 

naturally vanishing as ^ goes to zero. For simplicity 
we only show the conservative limit, requiring that no ef¬ 
fects are generated at observable scales in the CMB. On 
the other hand, due to the constraints from the Bullet 
Cluster discussed above, values of ^ much larger than 0.3 
are excluded for the complete range of q;d shown in the 
figure. 

In the case of atomic DM, there is an additional con¬ 
straint from the acoustic oscillation peak that results 
from the recombination of dark atoms [37]. This turns 
out to be very constraining on the parameter space of the 
atomic DM, but it still leaves some plausible parameter 
space for, e.g., the model proposed in [ 8 ]. 


Dark Halo Stability 


Tkin 


4g* \ -\/I35 m|j 


(18) 


where Ckin is the ratio of the dark and visible photon 
temperatures at the time of kinetic decoupling. The de¬ 
termination of the exact effect of the DM/dark radiation 
coupling on the CMB is beyond the scope of this let¬ 
ter. Here we will simply require that decoupling happens 
above Tkin > 640 eV, so that the DM/dark radiation cou¬ 
pling only affects multipoles above I > 2500 and is thus 
unconstrained by the Planck data. This will lead to a 
more conservative limit than what would be allowed by a 
more detailed analysis, but will be used here as a robust 
constraint. 

It should be noted that values near the lower end 
of this limit, or slightly below it, can help to allevi¬ 
ate the missing satellites problem [34, 35]. The cut-off 


The dark pair plasma will heat up and become viri- 
alized both in galactic DM halos and in galaxy clus¬ 
ter collisions. Similar to visible matter, it can dissipate 
heat through dark bremsstrahlung and through Comp¬ 
ton scattering off dark photons. The characteristic time 
scale for dark bremsstrahlung cooling for the plasma pa¬ 
rameters given above is [ 8 ] 

O 3/2„l/2 

Wn.. « 4 3 ^'' « 6.7 • 1019 yr, (19) 

16 npa^ 

where the dependence on the DM mass cancels for our 
model: cooling becomes more efficient for smaller DM 
mass, but this is offset by the smaller coupling constant 
required to produce the correct relic density. The char¬ 
acteristic time scale for dark Compton scattering is even 
larger and can be neglected for all reasonable values of 
the DM mass and coupling. 
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FIG. 1: The constraints on the parameter space of the dark 
pair plasma model. The black contours show the relic abun¬ 
dance of the dark plasma, as a fraction of the total DM abun¬ 
dance. The orange shaded region is disfavored by the kinetic 
decoupling constraint, and the lower limit of that region is fa¬ 
vored for alleviating the missing satellites problem. The blue 
shaded region is excluded, in absence of a UV completion, if 
we require that the Landau pole of the dark ?7(1) coupling is 
above the Planck scale. 


A thermally produced pair plasma can therefore not ef¬ 
ficiently dissipate heat, so that dark pair plasma halos do 
not collapse to disks. This is in contrast to atomic DM, 
consisting of asymmetric light and heavy dark fermions. 
Because of the non-thermal production mechanism, the 
dark coupling strength is here independent of the DM 
density. Atomic DM with light “dark electrons” can 
therefore cool down sufficiently fast to collapse to a dark 
disk within the lifetime of the Universe. However, the 
counter-streaming instabilities in the dark plasma might 
give rise to nontrivial substructure within the galactic 
DM halo. Conclusive treatment of this issue would re¬ 
quire numerical simulations, and is beyond the scope of 
this letter. 


DISCUSSION AND CONCLUSIONS 

In this letter we considered the possibility that a sub¬ 
dominant component of DM has long-range interactions 
and exists as a dark pair plasma. We found that a ther¬ 
mally produced dark pair plasma is necessarily collision¬ 
less, but will self-interact through collective plasma ef¬ 
fects. The possibility of forming collisionless shocks will 


then modify the dynamics of galaxy cluster collisions, 
leading to effects such as an offset of DM halos as in 
Abell 3827, or to an isolated DM clump as in Abell 520. 

To our knowledge, no numerical simulations or exper¬ 
iments have so far directly investigated non-relativistic 
collisionless pair plasmas. Our argument is however quite 
general, based on instability growth rates in the linear 
regime. For an atomic DM plasma which imitates the 
SM, the analogous behavior of visible astrophysical plas¬ 
mas can be directly used as a proof that collisionless 
shocks should also exist in the dark sector. Nevertheless, 
numerical studies of non-relativistic dark pair plasmas 
would be desirable for a more thorough treatment. 

Ultimately, collisionless shocks are an efficient form 
of DM self-interactions, explaining the features observed 
in the Abell 3827 and 520 clusters. If a galaxy moves 
through an intra-cluster medium of interacting DM, a 
drag force between the halo DM and the background clus¬ 
ter dark plasma is generated through plasma instabilities, 
potentially distorting halo shapes or even striping galax¬ 
ies of the interacting DM component. In cluster mergers 
bow-shaped dark matter shock fronts are expected to be 
observable in high resolution weak lensing studies. Thus 
a dark pair plasma offers spectacular signatures for its 
discovery, and we encourage the astrophysics community 
to look for these effects. 
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